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For a quantum-mechanical counting process we show er- 
godicity, under the condition that the underlying open quan- 
tum system approaches equiUbrium in the time mean. This 
implies equality of time average and ensemble average for cor- 
relation functions of the detection current to all orders and 
with probability 1. 

PACS numbers: 02.50.Ga, 02.70.Lq, 03.65.Bz, 42.50.Lc 



I. INTRODUCTION 

Modern research on quantum- mech anica l counting pro- 
be it numerical simulations 



[|Ca 



ar| or cxpcrimen- 



cesses, 

tal investigations [MYK|, usuaUy starts from the tacit 
assumption that for the study of statistical properties 
of the counting records it does not make a difference 
whether a large number of experiments is performed or 
a single very long one. This assumption amounts to er- 
godicity of these records. In several recent discussions, 
e.g. |BESW| , p?aS|jPlK| , pr4pCJv^ , investigators have ad- 
dressed the question of it s validity. A partial result was 
obtained by Cresser [Ore], who proved ergodicity in the 
L^-sense and to first order in the detection current. In 
this paper we establish ergodicity in the full sense (Theo- 
rem 3), in particular to all orders in the detection current 
and with probability 1 (Theorem 4). Theorem 5 formu- 
lates ergodicity in terms of multi-time coincidences. 

For the description of detection records we em- 
ploy the rigorous formulation of Davies and Srinivas 



| Da-v|JSrD |, whic h has set the tone for later investigations 
|Ca7,Wi^,|GaZ|. 



II. COUNTING PROCESSES ACCORDING TO 
DAVIES AND SRINIVAS 

We consider an open quantum system under continuous 
observation by use of a finite number k of detectors. The 
state of the system is described by a density matrix p on a 
Hilbert space, obeying a Master equation p = Lp, where 
L is a generator of Lindblad form |Lin]. Normalisation is 
expressed by the relation 



trL(p) = for all p 



(2.1) 



A counting process connected to this quantum evolution 
is based on an unraveling of the generator 



L = Lo + J2Ji 



(2.2) 



which is interpreted as follows. The reaction of the de- 
tectors to the system consists of clicks at random times. 
The evolution p i-^ e*^°{p) denotes the change of the 
state of the system under the condition that during a 
time interval of length t no clicks are recorded. The op- 
erator p ^ Ji (p) on the state space describes the change 
of state conditioned on the occurrence of a click of de- 
tector i. For computational convenience we assume these 
operators to be bounded. So, if p describes the state of 
the system at time 0, and if, during the time interval 
[0,t], clicks are recorded at times ti,t2, ■ ■ ■ ,tn of detec- 
tors ii,Z2, . . . ,in respectively, and none more, then, up 
to normalisation, the state at time t is given by 

g(t-t„)Loj^^g(t,.-t„-i)Lo...g(t2-ti)Loj^^gtiLo(p) (2.3) 



The probability density /*((ii,ii), 



,in)) for these 



cHcks to occur is equal to the trace of ( |2.3| ). 

We imagine the experiment to continue indefinitely. 
The observation process will then produce an infinite de- 
tection record ((ti, ii), (t2,*2), (^3,^3), • • ■) , where we 
assume that < ti < t2 < < . . . , and lim„^oo tn = 00 
(i.e., the clicks do not accumulate). 

Let VI denote the space of all such detection records. 
By an event we mean some property of the record, which 
we identify with the set E d Vl oi all records with this 
property. The events decidable at or before time t > 
form a cr-algebra Et |Dav |. Together these cr-algebras 
generate the full ti-algebra E. Following Davies and Srini- 
vas we may now formulate the effect of observation on the 
quantum system as follows: If t is a positive time, E an 
event in Ej, and p denotes a state, then we define 

Mt{E){p) := 



n—0 — 1 in—^ 



JO JO 



g(t-t„)Lo J. g(t„-t„_i)Lo ^ ^ ^ e''*^~*^'^^° Ji^e*^^° {p) 

xdtidt2---dtn . (2.4) 

Here 1e denotes the indicator function of the event E 
and Mt{E) is the effect on the quantum system of the 
occurrence of i? G Et . Then 



1 



;(£;) ■.= ivMt{E){p) 



(2.5) 



is the probability of the occurrence of E given that the 
system starts in p. We extend the notation (2.5) also to 
density matrices p which are not normalised. The count- 
ing process as a whole is described by the family (Mf )t>o. 
The effect of the counting on the quantum system, when 
the outcome is ignored, is the time evolution 

Tt{p) Mt{n){p) . 

It follows from the Dyson series ( p.4[ ) with E = ^ that 
Tt i s ind eed the original time evolution e*^, in particular, 
by (2J), Tt preserves the trace. 



III. ERGODIC THEORY 

The time shift by t seconds is described by the map 
at on r2, which is given on a particular record lo — 
((ii,ii), (t2,«2), (^3,^3), ...) e f2 with tk < t < tk+i 
by at{uj) := {{tk+i ~ t,ik+i), {tk+2 - ife+2), ■ • ■)• The 
time shift of an event E towards the future is given by 

The crucial property of the counting process {Mt)f>a 
is the following. For all s,i > and all events £' e S^, 
F G 'St we have 



Ms+tiF n a^\E)) = MsiE) o Mt{F) 



(3.1) 



This Markov property was proved in |Dav|. Putting E ~ 
F — we recover the semigroup property Tg+t = o Tt 
of the time evolution. 

When F e Et and s > then P^+'*(F) does not depend 
on s. Indeed, since n = (Jt^{n) and Ts preserves the 
trace, 

Fl+^iF) = tr {Mt+s{F)ip)) ^tv{Mt+siFna^\mp)) 
^ = ^ tr {Ms{n) o MtiF){p)) = tr (T, o Mt(F)(p)) 
= tr{MtiF)ip))^VliF). 
Therefore, by Kolmogorov's extension theorem, the fam- 



iiy 



t>o of probability measures on the a-algebras 



(St)t>o with densities (/*)t>o extends to a single proba- 
bility measure Pp on the full cr-algebra S. 

Lemma 1. For all t > 0, all E e T,, F e T,t, and all 
states p: 

Pp{Fn<j^\E))^FM.iF)ip){E) . (3.2) 
In particular, 

Fp{<j^\E))^FtUE) ■ (3.3) 

Therefore, if p is invariant under Tt, then Pp is a station- 
ary probability measure on fl. 



Proof. First suppose that £' e E^. Eq uality (3.2) is 
obtained from the Markov property ( b.lh b y acting on 
p and taking the trace on both sides ( p?.3D follows by 
putting F = n. The statements extend to all G E by 
Kolmogorov's extension theorem since s was arbitrary. □ 



Definition. 

• The evolution {Tt)t>o of a quantum system is said 
to converge in the mean to an equilibrium state p 
if for all normalised density matrices 'd and all ob- 
servables x: 



lim — 



tr ((Tti?)x) dt = tr (px) 



• The counting process {Mt)t>o will be called ergodic 
if the following holds. Given any time-invariant 
event E, i.e. cTt^{E) — E for all t > 0, then either 
Fi){E) = for aU density matrices d or F^{E) = 1 
for aU 

The condition on (rt)f>o is satisfied in many cases of 
practical importance. 

Theorem 2. If the evolution Tt = e*^, t > 0, converges 
in the mean, then t he co unting process {Mt)t>o is ergodic 
for any unraveling (2.2 ). 



Proof. Le t E be a time- invariant event and i9 any state. 
Then by (fj), P{,{E) = ¥4a:^\E)) = FtME)- Since 
P^ is linear and continuous in i?, we may average both 
sides over the interval [0, r] and take the limit t —>■ 00 
to obtain P,g{E) = Pp{E). For an unnormalised density 
matrix x we find instead that 

P^{E)=P,{E)tr{x). (3.4) 
If F is any event in Ef then 

^^{FnE) ^ ¥4Fna^\E)) * = ^ PM,(F)w(i^) 
Pp(£;) tr {Mt{F){^)) = Fp{E)Fi){F) 
F^>{E)F4F) . 



(3.4) 



(3.4) 



The resulting equation extends to all G E, in particular 
it holds for F = E: 

F^iE)=F4Ef . 
It follows that F§{E) is equal to or 1. □ 

Let us denote the expectation f{uj)dFp{uj) of an in- 
tegrable function / on il by Ep(/). 

Theorem 3. If the evolution {Tt)t>Q converges in the 
mean to p, then for all integrable functions h on and 
all initial states d we have, almost surely with respect to 



lim 



h{at{uj))dt = Ep{h) 



(3.5) 



Proof. By Lemma 1 and Theorem 2, Pp is stationary and 
ergodic. Hence, by Birkhoff 's individual ergodic theorem, 
the limit on the left exists almost surely with respect to 
Pp, and is equal to the constant Ep(/i). Since the set F 



of points w G for which (3.5) holds, is time-i nvar iant , 
we have Fs{F) Fp{F) = 1 for all states ^ by (fj). □ 
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IV. APPLICATIONS 

The main result of the present ergodic theory for quan- 
tum counting processes, Theorem 3, can be made consid- 
erably more concrete by applying it to detection currents 
and multi-time coincidences, showing bunching or anti- 
bunching. 

For simplicity we consider only one detector, which re- 
sponds to a point event at time s by producing a current 
j{t — s) at time t. (This will be zero for t < s.) The total 
detection current is given by 



Let Pp be the unique stationary extension of Pp to 
negative times on the configuration space Q of the full 
real line. We shall denote expectation with respect to 
this measure by Ep. 

Theorem 4. Let the quantum evolution (Tt)t>o con- 
verge in the mean to a state p and let the detector re- 
sponse function 7 : M — ^ [0, 00) be bounded and inte- 
grable. Then for aU < ti < t2 < ■ ■ ■ < tn and all initial 
states 1) we have, almost surely with respect to F^, 

lim - / ItM^) ■ ■ ■ Iu+t{oj) dt = Ep {It, ■■■It J . 

For n — 2 this theorem implies a quantum-mechanical 
version of the Wiener-Khinchin theorem. In the proof we 
shall make use of the n on-exclusive pr obability density of 
the stationary process | vKa , paZ , Cre |, 

gniti,t2, . . . ,t„) := tr {jTt„^t„^iJ ■ ■ ■ JTt^-tiJ{p)) ■ 

The functions are related to the probability density 
/* from (2^) of the counting process (where i > tn), by 



gn{ti, t2 
ft 



m— 1 



/m+n({^l 



,tm} U {Si, 



J) 



xdsi ■ ■ ■ dsj, 



fi{ti,t2,...,tn}UUj)dLJ; (4.1) 



here fit is the set of finite subsets of [0,t], which can 
be identified with the time-ordered points in {0} U 
Um=i I*-*' ^]™- mean dsids2 ■ ■ ■ ds„i if a> = 

{Sl, S2, • ■ • , Srn} with Si < S2 < • • • < S,„. 

Proof of Theorem 4. First we note that Theorem 3 also 
holds if rj, Pp and Ep are replaced by i7, Pp and Ep re- 
spectively, as introduced above, and at by the left shift 
of w C M. Then we have /s+t(ijj) = Is{at{co)). Now fix 
n 6 N and < ii < • • • < i„. Let /i : f2 — > R be given by 

h{co) ■.^It,iL0)It,iLu)^^^ItjL0) . 

It follows that ho ut — It-,+tIt2+t ■ ■ ■ It„+t, and the state- 
ment to be proved follows from Theorem 3, provided that 
h is integrable. In the Appendix we shall show that this 
is indeed the case. □ 



As our second application we shall show that the non- 
exclusive probability densities gn have a straightforward 
pathwise interpretation: they are equal to the frequency 
of multi-time coincidences on almost every detection 
record. For this, let N[a^b]{^) ■— #(^ l~l 1'^^^]) denote 
the number of clicks detected during the time interval 
[a,b]. 

Theorem 5. Let {Tt)t>o converge in the mean to the 
equilibrium state p. Then for all n € N, all < ti < t2 < 
. . . < tn, all e between and mini<j<„(tj+i — tj), and 
all initial states we have, almost surely with respect to 



1 /■'^ / " 
lim — 

r — )-oo T 



g{si, . . .,Sn)dsi ■ ■■dsn . (4.2) 



Proof. Fix 71 G N and a sequence < ti < t2 < ■ ■ ■ < tn 
of times. Let : O — > {0, 1} be the function that maps 
cij e r2 to 1 if w contains exactly n points, one in each of 
the intervals [ti, ii+e], • • • , [tn, tn-\-e], and to otherwise. 
Then we obtain for t > tn + e, u sing set notation and the 
integral-sum lemma from [ LiM |, 



tl+S 



, Sn) dsi • • • dSn 

(4.1) 



K{a)g{a)da"^^' j / K{a)f\aiJ (3)dad(3 



(4.3) 




K{a) f\io)dw. 
' I 

A short calculation shows that 

n 

^if(a) = n%.^*.+-lH- 



(4.4) 



Since < 5„(si,S2, . . . ,s„) < the integral (]4.3| ) is 

convergent, hence the product on the r.h.s. of (4.4) is 
integrable as a function of uj. Application of Theorem 3 
to this product now yields the statement. □ 



V. DISCRETE TIME 

There is an obvious anal ogue o f our main result (The- 
orem 3) in discrete time [ MaK | . A Kraus measurement 
[ Kra I is given by a decomposition of a completely positive 



operator T on state space as 



3 



where p ^ aipa* describes the state change of the den- 
sity matrix p when the measurement gives the outcome i. 
Thus for initial state i? the probabihty of finding the se- 
quence of outcomes zi, i2, . . . , ?m by repeated Kraus mea- 
surement is given by 

tr (ai„ • • -a^^da*^ ■ ■ ■ a*^) . 

As in continuous time, this yields a probabihty measure 
Vi) on the space of detection records J7 := {1, 2, • ■ • , k}^ . 
Again, if (r")„gN converges in the mean to some state 
p, then the only time invariant events in SI have measure 
or 1 for all P^. In particular, Pp is ergodic. 



APPENDIX: 

We shall show that, in the situation of Theorem A, h :— 
Iti ■ ■ ■ /*„ is an integrable function on provided that the 
jump operator J is bounded and the detector response 
function 7 : M — > [0, 00) is bounded and integrable. 

Let M max(l, ||7||oo)- Fix ?i e N and a sequence 
< ti < <2 < • ■ • < in of times. Let 



00 

m=0 



(p{si - t) • • • ip{Sra - t)dsi ■■ - dSn 



Therefore, since the r.h.s. does not depend on 

Ep(/ti ••■/*„)- lim IEp((/t, ■■■Itjoat) < 00. 



□ 
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'Pit) ■■=Y.^{tj-t) ■ 

Then ip is also integrable, with \\ip\\i = rt||7||i. For G N, 
let Jn.k denote the set of all surjections {1, •■•,«} 
{1, ■ ■ ■ , k}. Then we may write for any w G f2, 

= ^ ■ ■ ■ ^ litl - Si)- ■ ■■^{tn- Sn) 

n 

= E E E 7(ii-aj(i))---7(in-aj(n)) 



k=l 



< n ■ n' 



(Al) 



Using set notation and the integral-sum lemma |LiM] 
again we conclude that, for all t > and u > tn + 1, 

Ep((/tjt, •••/tJoat)/Af"n"+i 

(Al) 



•^f^" aCuj \sea ) 



f2„ JQ. 
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